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Abstract 



The concept of the theory of continuous groups of transformations has attracted the 
attention of applied mathematicians and engineers to solve many physical problems 
in the engineering sciences. Three applications are presented in this paper. The first 
one is the problem of time-dependent vertical temperature distribution in a stagnant 
lake. Two cases have been considered for the forms of the water parameters, namely 
water density and thermal conductivity. The second application is the unsteady free- 
convective boundary-layer flow on a non-isothermal vertical flat plate. The third 
application is the study of the dispersion of gaseous pollutants in the presence of a 
temperature inversion. The results are found in closed form and the effect of param- 
eters are discussed. 

1 Introduction 

According to Seshadri and Na [1], different methods for carrying out similarity analysis of 
partial differential equations are classified into direct methods, where the concept of group 
invariance is not explicitly invoked, and group-theoretical methods that are based upon 
the invocation of invariance under groups of transformations of the partial differential 
equation and the auxiliary conditions. 

The foundation of the group-theoretical method is contained in the general theory of 
continuous transformation groups that were introduced and treated extensively by Lie [2] 
in 1875. Group-theoretical methods provide a powerful tool because they are not based 
on linear operators, superposition, or any other aspects of linear solution techniques. 
Therefore, these methods are applicable to nonlinear differential models. 

Throughout the history of similarity analysis, a variety of problems in science and 
engineering have been solved. Among these we find a general treatment of steady two- 
dimensional incompressible laminar flow of fluid into an infinite region of the same fluid 
by Abbott and Kline [3] in 1960, impact of thin long rods by Taulbee et al [4] in 1971, 
wave propagation in viscoelastic, viscoplastic and electrical transmission lines by Ames 
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and Suliciu [5] in 1982, time-dependent free surface flows under gravity by Sachdev and 
Philip [6] in 1986, problem of ocean acoustics by Richards [7] in 1987, unsteady free- 
convective boundary-layer flow on a non-isothermal vertical flat plate by Abd-el-Malek 
et al [8] in 1990, steady free-convective boundary-layer flow on a non-isothermal vertical 
circular cylinder by Abd-el-Malek and Badran [9] in 1991, dispersion of gaseous pollutants 
in the presence of a temperature inversion by Badran and Abd-el-Malek [10] in 1993, and 
nonlinear temperature variation across the lake depth neglecting the effect of external 
heat sources by Abd-el-Malek [11] in 1997. Recently Boutros et al [12, 13] considered two 
problems namely, potential equation in triangular regions with temperature distribution 
along the boundaries in the form of polynomials, and the second is the time-dependent 
vertical temperature distribution in stagnant lake taking into consideration an external 
heat source. Many physical application are illustrated in Sedov [14], and Rogers and 
Ames [15]. 

In the present paper we will present three physical applications of the group-theoretical 
method namely, time-dependent vertical temperature distribution in a stagnant lake tak- 
ing into consideration the effect of external heat source. The problem has been solved for 
two possible forms of the water parameters (water density and thermal conductivity). Sec- 
ond application is the unsteady free-convective boundary-layer flow on a non-isothermal 
vertical flat plate. As a third application we study the dispersion of gaseous pollutants 
in the presence of a temperature inversion. The obtained results are found in closed form 
and the effect of different parameters are discussed. 

2 Application (I): Time-dependent vertical temperature 
distribution in a stagnant lake 

2.1 Mathematical formulation 

Consider the one dimensional heat transfer equation for heat flux in the vertical direction, 
neglecting the convective motion of the fluid and assuming the absolute value of the 
specific heat of the water as sensibly constant within the range of temperatures considered; 
take it to be unity. The vertical transport of heat in a deep lake is modelled by 



where "T" is the temperature; r(z,t) is rate at which solar radiation is absorbed by the 
water; "t" is the time; "z" is the distance measured downward from the lake surface; "p" 
is the density; and "re" is the thermal conductivity. 

Boundary and initial conditions: 

During early Spring, most of the lakes exhibit a nearly homothermal temperature distribu- 
tion with a low degree of temperature extending all the way to the bottom (see Sundaram 
and Rehm [16]). In all of the calculations presented in this paper, the initial condition 
will be taken as that corresponding to the end of spring homothermy, i.e., 



p(T)T t = [K(T)T z ] z + r{z,t), 



z>0, t>0, 



(2.1) 



T(z,0) =T 



(2.2.1) 



where Tq is the temperature of the lake at maximum Spring homothermy. 
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The boundary conditions are considered as follows: 

(i) T z (0,t)=0, t>0, (2.2.2) 

(ii) T(h,t) = T + 7t 1/m , t > 0, < 7< 1, m > 0. (2.2.3) 

In our analysis we will consider two cases for the form of the density p(T) and the 
thermal conductivity k(T), namely: 

Case (1): p = aq{z)(T - T ) m , K = {3g(z), 

Case (2): p = aq{z)(T - T ) s , K = (3(T-T ) n , 

where a and /3 are constants, m, s and n are positive constants and q(z) and g(z) are 
arbitrary functions to be determined later on. 
Write 

T(z,t) = w{z,t)+T , (2.3) 
by which differential equation ( |2.lD takes the form: 

p(w)w t = [k(w)w z ] z + r(z,t), z>0,t>0, (2.4) 

and the initial and boundary conditions take the form: 

(1) initial condition: 

w(z,0)=0, (2.5.1) 

(2) boundary conditions: 

(i) w z (0,t)=0, t>0, (2.5.2) 

(ii) w(h,t) ='jt l/m , t>0, 0<7<1, m>0. (2.5.3) 

2.2 Case (1): p = aq(z)w m , k = j3g{z) 

In this case differential equation ( |2.4| ) takes the form: 

f3g(z)w zz + (3g z w z - aq{z)w m w t = -r(z, t). (2.6) 

2.2.1 Solution of the problem 

Our method of solution depends on the application of a one-parameter group transforma- 
tion to the partial differential equation ( |2.6| ). 

2.2.1 (a) The group systematic formulation. The procedure is initiated with the 
group G; a class of transformation of one-parameter "a" of the form 

G : Q = C Q (a)Q + P Q (a), (2.7) 

where Q stands for t, z, r, w, k, p and the C's and P's are real-valued and at least 
differentiable in "a". 
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2.2.1 (b) The invariance analysis. To transform the differential equation, transfor- 
mations of the derivatives of w, k and p are obtained from G via chain-rule operations: 



Oil 



Qij'i 



z,t; j = z,t. 



(2i 



Equation ( |2,6| ) is said to be invariantly transformed, for some function H{a), whenever 



09 Wzz + Pg z w z - aq(w) m w T + r = H(a) [0gw zz + f3g z w z - aqw m w t + r] 
Substitution from ( |2.7| ) into ( |2,9| ) yields 







C g C v 



gwzz + P 



C g C v 



(C 



z\2 



g z w z - aw 



qwt + [C r ]r + R(a) 



(2.9) 



(2.10) 



H{a) [(3gw zz + (3g z w z - aqw m w t + r] 



where 



R(a) 



ft 



pg C v 



(C 



z\2 



aP q {C w w + P w ) m ^— 



k=l 



m 
k 



<x{C q q)—w t Yj( k I (C W w) m - k (P w ) k + P r . 



The invariance of ( ^S| ) implies R(a) = 0. This is satisfied by putting 
pr _ pq _ pw = pg = 

and 

C r = H{a). 



~ (jq^jw^m+l ' 




- C g C w- 









(2.11) 



(2.12) 



Moreover, the boundary and initial conditions ( 2.5.2|) , ( |2.5.3 ) and (2.5.1) are also in- 
variant in form, implying that 



pz 



0, P T = 0, C z = 1 and (C l 



Combining equations (2.12) and invoking the result (|2.13j ), we get 
C r = C q C w = C 9 C W which yields C q = C 9 . 



(2.13) 
(2.14) 



Finally, we get the one-parameter group G which transforms invariantly the differential 
equation fl2.1|), as well as the boundary and initial conditions (2.2). The group G is of the 
form 

G : z = z, t={C w ) m t, q = C q q, r = C q C w r, w = C w w, g = C q g. (2.15) 

2.2.1 (c) The complete set of absolute invariants. Our aim is to make use of group 
methods to represent the problem in the form of an ordinary differential equation. Then 
we have to proceed in our analysis to obtain a complete set of absolute invariants. 
If r/ = rj(z, t) is the absolute invariant of the independent variables, then 



gj(z,t;w,r,K,p) = Fj[rj(z,t)]; j = 1,2,3,4 



(2.16) 
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are the four absolute invariants corresponding to w, r, k and p. The application of a basic 
theorem in group theory, see [17], states that: a function g(z,t;w,r, k, p) is an absolute 
invariant of a one-parameter group if it satisfies the following first-order linear differential 
equation 



i=l 



dg_ 

dQi 



0; 



Qi = z,t,w,r,n,p, 



(2.17) 



where 



dCQ \ 0, a 



QpQi 

da 



1,2, ... ,6, 



(2.18) 



and a denotes the value of "a" which yields the identity element of the group. 
Owing to equation (|2 . 1 7|) , r](z,t) is an absolute invariant if it satisfies 

(aiz + (5i)n z + (a 2 t + 2 )r]t = 0, 
Group ( 2.15| ) gives: 

«i = A = #2 = 0, 
and hence from (2.17) and ( 2.2C| ) we get 

Vt = 0, 
which gives 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



r](z,t) =F(z). 

Without loss of generality we can use the identity function: 

7]{z,t) = z. (2.23) 
By a similar analysis the absolute invariants of the dependent variables w, r, k and p are 



w(z,t)=T(t)F(r ] ), r(z,t)=A(t)Q( V ), 
q(z) = B(t)$( V ), g(z) = Y(t)9( V ). 

From which we conclude that 
q(z) = $(??), 

g(z) = ^(r]). 



At t = : T(0) = 0. 



(2.24) 

(2.25) 
(2.26) 
(2.27) 
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2.2.2 The reduction to ordinary differential equation 

Substituting from ( p. 24 ) into ( |2.6| ) and dividing by T, we get 

0VF m + P9 V F V - a<S>T m - l F m+l T t = _ A ^^\ ( 2 .28) 

For ( p. 28 ) to be reduced to an expression in the single independent invariant 77, it is 
necessary that the coefficients should be constants or functions of 77 alone. Thus 

r™- 1 ^ = Ci, 4^ = C 2 . (2.29) 
Take C x = 1: 

r(t) = (mi) 1/m , (2.30) 
(gjH ) and (|2~30l) yield: 

A(t) = C 2 (mi) 1/m , (2.31) 
Hence fl2.28| ) may be rewritten as 

/3*F W + f3^ n F v - a<5>F m+1 = -C 2 @(rj). (2.32) 
Following Girgis and Smith [18], we assume the heat source distribution in the form: 

8(77) = e"^, (2.33) 

where £ is the absorption coefficient for water, which has the value 0.048. 
Take: 

® (v) = Wj> F (^ ' 0<il<h. (2.34) 
According to ( |2.33 ) and ( 2.34j ), equation (|2.32j ) takes the form 

p-VF m + p^E,, - a*F m = -C 2 e-^. (2.35) 
Write: 

= e~ w , (2.36) 
where \i is a constant, then (|2.35| ) becomes 

F m - ^ - ^F m = -^ e -«-">", (2.37) 

with the boundary conditions: 

(i) F(0) = 0, (2.38) 

(ii) F(h) =j/m (1 / m l (2.39) 
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2.2.3 Analytical solution for different forms of the parameters 



For m = 1, differential equation (2.37) becomes 



G 







cr 



a 




and the boundary conditions become 

(i) F(0) = 0, 

(ii) F(h)= 7 , 

which has the exact solution 

F( v ) = a ie riri + a 2 e r2V + a 3 e~^-^ , 

where 



(2.40) 

(2.41) 
(2.42) 

(2.43) 



r l,2 



fi ± yV + 4(J 2 



ri for (+) sign and r 2 for (— ) sign, and 
Ch 

>( M -e) 2 -/i/3(/x-6-«' 



«3 



For finite temperature, a\ = 0, and applying condition ( |2.41| ), we get the solution in 
the form: 



F{rf) = a 3 



Q-iS-riv _|_ * £ e W/ 

7*2 



Hence the temperature distribution across the lake, corresponding to case (1) is: 



T(z,t)=T 



Cot 



/3fr - £)2 _ ppfr _ £) _ Q 



^2 



(2.44) 



(2.45) 



For < i < 150 (in days), following Girgis and Smith [18], we use the following values 
of the parameters: C 2 = 2496, /? = 12355, p = 1.439239 x 10~ 4 , £ = 0.048, h = 400 meter, 
Tq = 4°C. The obtained results are plotted in Fig.l and Fig. 2. 



2.3 Case (2): p = aq(z)w s , k = (3w n 
Differential equation fl2.4| ) takes the form 

w n w zz + nw n ~ 1 {w z ) 2 — a 2 q(z)w s Wt = —r(z, t). 

Following the same analysis as in case (1), we get the following group G: 

G: z = z, t = {C w ) m t, q = {C w ) n+m - s q, r = (C w ) n+1 r, w = C w w. 



(2.46) 



(2.47) 
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Figure 1. Distribution of temper- 
ature T (time = 40 days) againts 
the lake depth "z" in meters, cor- 
responding to case: p = aq(z)w, 
k = (3g{z) for different values of pa- 
rameter "a". 
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Figure 2. Distribution of tempera- 
ture for different times to constant "a" 
(a = 14095), for p = aq(z)w, k = 
f3g(z). 



The absolute invariants are: 

r, = z , w(z,t)=T{t)F(ri), r(z,t) = A(t)Q{r]), q(z) = B(t)^(rj). 
Again, it is clear that B(t) = 1; from which we get 

q(z) = $(77), 
leading to 

Q = q, 



(2.48) 



(2.49) 



(2.50) 



which is satisfied if and only if (c w ) n + m ~ s = 1. That is C w = 1 or n = s - m. C w can 
not be unity. Hence we obtain the only possible case: 



n = s — m. 



Take 



$(77) 



1 



F(r,Y 

At t = : T(0) = 0. 



F(ri) + 0, < r) < h, 



(2.51) 

(2.52) 
(2.53) 



322 



M.B. Abd-el-Malek 



4.00 



0.00 



40.00 



80.00 



120.00 " 




0.00 



40.00 



80.00 



120.00 



5.00 
I 



15.00 25.00 35.00 °C 

J I I I 







' / 
' / 
■ '/ 
'/ 


40 days 




120 days 




150 days 



Figure 3. Distribution of temper- 
ature T (time = 40 days) againts 
the lake depth "z" in meters, cor- 
responding to case: p = aq(z)w, 
k = (3 for different values of param- 
eter "a". 



Figure 4. Distribution of tempera- 
ture for different times to constant "a" 
(a = 13306), for p = aq(z)w, k = 
(3g(z). 



Following the same analysis as in section (2.2.2), we reach to the following ordinary 
differential equation: 

F n F m + uF"" 1 ^) 2 - a 2 F s = -C 2 e~^. (2.54) 
For the case when m = s = 1, and from ( |2.51| ), we find n = 0. Hence (2.54) becomes 



wo 



a 



! F = -C 2 e"^, 



(2.55) 

(2.56) 
(2.57) 

a 2_ e - ■ ( 2 - 58 ) 

For finite temperature, a\ = 0, and applying boundary condition fl2.56|) , we get the solution 
in the form: 



and the corresponding boundary conditions are: 

(i) F(0) = 0, 

(ii) F(h)=r 

Applying boundary conditions Q2.56|) and ( |2.57D , we get the solution 

F{ v ) = aie^ + a 2 z- un + 



Co 



a 2 -e 



-CTT) 



Hence the temperature distribution across the lake, corresponding to case (2) is: 



T(z,t) =T + 



C 2 t 



a 



(2.59) 



(2.60) 
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For < t < 150 (in days), following Girgis and Smith [18], we use the following values 
of the parameters: C 2 = 0.2014, = 12391, £ = 0.048, h = 400 meter, T = 4°C, the 
obtained results are plotted in Fig.3 and Fig.4. 



3 Application (II): Unsteady free-convective boundary- layer 
flow on a non- isothermal vertical flat plate 

3.1 Mathematical formulation 

Consider a natural convective, laminar, boundary layer adjacent to a semi-infinite, vertical 
flat plate. The fluid is isothermal and of constant temperature T^, far from the plate. The 
plate has nonuniform surface temperature T w > (heated plate case). The fluid has 
the following constant properties: "/?" is the volumetric coefficient of thermal expansion, 
V the kinematic viscosity, and "a" is the thermal diffusivity. 

Along with the application of the Boussinesq and boundary-layer approximation, the 
equations of motion may be written as: 

(1) Conservation of mass: 

Ux + v w = 0. (3.1) 

(2) Momentum equation: 

u I + uUx + vUy = gl3{T-T 00 ) + vUyy. (3.2) 

(3) Energy equation: 

T- t + uTx + vT v = aTyy. (3.3) 
Boundary conditions: 

(i) U = 0, u = 0, T = T w (x,t) at y = 0, t > 0, (3.4.1) 

(ii) u = 0, T = Too as t/^oo. (3.4.2) 
Dimensionalize the variables according to: 



w 



u = -, v=(Gt)/-, t = U v 



(3.5) 



where L is some arbitrary reference length, AT = T re f — T^, T re f is some arbitrary 

reference temperature, U = (gPLAT) 1 / 2 is the Characteristic velocity, and Gr = 

is the Grashof number. 

In dimensionalized form 

The basic equations are: 



U X +Vy = 0, 



(3.6) 
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U t + UU X + VUy = T + Uyy, (3.7) 



T t + uT x + vT v = —T vv ; Pr = - is the Prandtl number. (3.8) 
Pr yy a 



The boundary conditions are: 

(i) v = 0, u = 0, T = T w {x,t) at y = 0, t > 0, (3.9.1) 

(ii) u = 0, T = 0, as y -► oo. (3.9.2) 
From the continuity equation there exists a stream function ty(x,y) such that 

— —■ fc- (3 ' 10) 

Momentum and energy equations take the form: 

Vyt + ^y^yx ~ ^x^yy = @T W + ty yyy , (3.11) 

T w o t + e(T w ) t + T w * y e x + e^(r w ) x - r^e^ = ^T w e yy . (3.12) 

Boundary conditions are: 

V x (x, 0, t) = V y (x, 0, t) = 0, 0(x, 0, t) = 1, 



lim ty v (x,y,t) = 0, lim @(x,y,t) = 0. 

y— >oo y—*oo 



(3.13) 



3.2 Solution of the problem 

A class of two-parameter group (01,02) has the form 

S = C S (a 1 ,a 2 )S + K s (a 1 ,a 2 ), (3.14) 

where U S" stands for x, y, t; tp, 0, T w , and C's and K's are real valued and differentiable 
functions with respect to a\ and a 2 . 

The invariance analysis: 

Invariance of the transformed momentum equation leads to the following group G: 



G : I 



x = (C y C^)x + K x , y = C y y, t = {C y ft + K l 



v = c*v + k*, t w = [-^—)t w , e = e. 



The absolute invariants: 

A function rj(x, y, t) is an absolute invariant of a two-parameter group if 77 satisfies: 

(ctix + a 2 )r] x + (a 3 y + 0:4)% + ("5* + a 6 )r}t = 0, and 

(3.16) 

(fax + p 2 )vx + {fay + Ai)% + (fet + A>H = 0, 
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where 

. Q Q. 8K X . Q Q. 

a x = — — {<q,a%), ol 2 = — — [al,a%), 
oa\ oai 

ft = ^(<M), & = ^(a?,a°), .... 
aa2 oa<i 

The only possible form for the absolute invariant rj is: 
y/ait + 6i ' 

where ai = as = &i = «6 = 06 are constants. 



Abd-el-Malek, ei al [8] have shown that the case of r] = y/\/a\x + bit + c does not lead 
to any solution. 

The complete set of absolute invariants: 

The only possible forms for the absolute invariant \& and T w are: 

*(x,y,t)=T(x,t)F{r,), (3.18) 
T w =T u(x,t). (3.19) 

3.2.1 Solution corresponding to the form of 77 in ( |3.17| ): 

The corresponding differential equations are: 

F vvv + (y»7 + F ) Fm ~ ^ + aiFr > + = °' 

1 (3 - 20) 

p^e w + (y 77 + Fj Q v + (2a x - F V )Q = 

with the boundary conditions: 

F(0) = F v (0) = 0, 6(0) = 1, F v (oo) = 0, 8(oo) = 0. (3.21) 
We get the following solution: 



x + 6 2 x + 62 „ 



(0.4472ait + 6i) 2 ' att + W 



(3.22) 



+ z/ (ait + 6 2 ) 5 /2' 



3.3 Conclusion 

It is clear from the obtained results ( |3.22j ) that the temperature profile overshoots in the 
region of the boundary layer near the plate. This phenomena occurs for values of a\ > 1, 
and becomes stronger as a% increases. This means that this phenomenon is accomplished 
by those cases for which T w decreases rapidly with time. 

If we study the effect of Pr on the temperature profile we find that there is a rapid 
increase in 9 near the plate. This becomes more evident for larger values of Pr. Also the 
thermal boundary-layer thickness decreases for increasing values of Pr. 
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4 Application (III): Dispersion of gaseous pollutants in the 
presence of a temperature inversion 

4.1 Mathematical formulation 

The gaseous pollutant is bounded from above by the ground surface and from below by 
the inversion layer, which is at height u h n from the ground surface. Assuming that the 
pollution, with concentration C(x,y), is evenly distributed throughout the layer, and the 
mean concentration of the pollutant at x = averaged over < y < h is constant and 
equal to Co . The diffusion of the pollutants takes place due to the wind that has a constant 
mean velocity u = u{x) in the x-direction, and the eddy diffusivities k\ and k 2 in the x 
and y-directions, respectively, are also independent of y. 

The normalized steady state diffusion equation, that governs the dispersion of the 
gaseous pollutants is 

u(x)C x = Kl(x)C xx + K 2 (x)C yy , (4.1) 

with the boundary conditions 
K\Cy = X-fC at y = 

C y = at y = 1 I 2 h 




7 = — , (4-2) 
C = 1 at x = 0, 0< y < 1 | u o 

c x = o 

where all x and y are scaled with respect to h, C with respect to Co, u with respect to no, 
Ki and K2 with respect to uoh, and no is a reference velocity. Values of A classify two cases, 
case (1): A C 1 corresponds to the case where no pollutant is absorbed by the ground, 
case (2): A > 1 corresponds to the case where all pollutant is absorbed by the ground. 
Introduce the non-dimensional function 9(x,y) and C*(x) such that 

C(x,y) = e(x,y)C*(x), (4.3) 

equation fl4.1| ) becomes 

u(c*e x + ec*) = K X {c*e xx + 26 x c* + ec* x ) + K 2 c*e yy . (4.4) 

4.2 Solution of the problem 

Following the same analysis as we did in application (I), we find that the group G\ which 



transforms invariantly the differential equation (|4.4| ) and the boundary conditions ( |4.2[) , 
is in the form: 

(x = E x (a)x, y = y, u = E x (a)u, 

1 ' {K 1 = (E x (a)) 2 K 1 , k 2 = k 2 , C* = E c *(a)C*, 6 = 9. 

The absolute invariants of the independent and dependent variables are: 



r] = y, F(x, u) = -, G(x, Ki) = -i, 
x x z 

= 9(y), k 2 = k 2 (x), C* = C*(x). 



(4.6) 
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4.3 The reduction to ordinary differential equation 

Substituting from ( |4.6| ) into ( |4,4| ) gives: 

[K2C*}8 yy + [x 2 GC* x - xFC*}8 = 0. (4.7) 

The requirement of reducing (|4.7j ) to two ordinary differential equations, for some constant 
p 2 , implies that: 

[x 2 GC* xx -xFC* x ]=p 2 [ K2 C*}. (4.8) 

Under this assumption, ( [4.4| ) gives the ordinary differential equation of 9(y), namely 6 yy + 
p 2 6 = 0, which has the solution 

0(y) = Acosp(y + e). (4.9) 
Rearranging fli~8|) , we get 

aC* xx -2pC* x -p 2 C* = 0, (4.10) 
where, using (|4.6| ), the constants a and /3 are 

x2 ° Kl a a xF u /a 

a = = — and p = - — = - — . (4-11) 

The ordinary differential equation ( |4.10| ) has the general solution 

C*(x) = A ie mix + A 2 e m2X . (4.12) 
Applying the boundary conditions ( |4.2j ) as x — > oo, for a > 0, we get the solution 

C*(x) = A ie mx , (4.13) 

where 



m 



P ~ V/3 2 + ap 2 . (4.14) 

Substituting (|]|) and (|1|) in (U), we get 

C(x,y) = Be mx cos p(y + e). (4.15) 

Application of the boundary condition at the inversion level determines e = — 1. 

Hence 

C(x,y) = Be mx cos p(y-l). (4.16) 
Application of the boundary condition ( (4.2| ) at the ground surface yields 

tanp=— . (4.17) 

K2P 

The constants "5" and "p" will be determined corresponding to the two limiting cases 
of A. 
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4.4 Analytical solution corresponding to the two limiting cases of A 

Case (1): A <C 1, which corresponds to no pollutant absorbed by the ground. 

A7 

Assuming to be a very small quantity leads to 

p 2 = * (4.18) 

and hence, by applying the boundary condition (|4.2j), we find B = 1. 
The concentration distribution for the case (1) is 

C = e mx . (4.19) 



where "m" is given by ( 4.14 ). 

Case (2): A> 1, which corresponds to all pollutant absorbed by the ground. 
In this case, from equation ( 4.17 ), we get 

p = Ntt. 

The concentration distribution model will be 



C(x,y) =Y,B n e mx cas[Nir(y-l)], 



71=1 



where 



m 



N 



2n- 1 



n = 1, 2, 3, . . . . 



(4.20) 



(4.21) 



The constants B n will be determined as the cosine's Fourier coefficients of the expansion 
of the function: C=lin0<y<l and x = 0, we get 



B„ 



sin Nir 
Ntt 



Hence the concentration distribution for the case (2) is 



71=1 



sin Ntt 
Ntt 



e mx cos[7V7r(y-l)], 



(4.22) 



where "m" and "iV" are given by (4.21). 
4.5 Results and discussion 

For the case where no pollutant is absorbed by the ground surface, it is found that the 
concentration distribution has the form 



C = exp 



l-L 



u 



■\/ u 2 + 4A/?i7 



which is independent of "y" and mainly depends on u, k\, k%, h and A. For very small A 
or very large u, it is clear that no pollutant will be absorbed by the ground surface. Also 
we concluded that: 
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(1) As "u" increases, the absorbed pollutant by the ground surface will be less. 

(2) As the eddy diffusivity ratio, a, increases the absorbed pollutant by the ground 
surface will be more. 



For the case where all pollutant is absorbed by the ground surface, it is found that the 
concentration distribution has the form: 



C 



2 E 

71=1 



sin Ntt 
Att 



cxp 



— (u - \/u 2 + 4N 2 TT 2 Kl) 

2k2 V " / 



cos[A7r(y — 1)]. 
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